Finite temperature phase diagram of a polarized Fermi gas in an optical lattice 
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We present phase diagrams for a polarized Fermi gas in an optical lattice as a function of tem- 
perature, polarization, and lattice filling factor. We consider the Fulde-Ferrel-Larkin-Ovchinnikov 
(FFLO), Sarma or breached pair (BP), and BCS phases, and the normal state and phase separation. 
We show that the FFLO phase appears in a considerable portion of the phase diagram. The dia- 
grams have two critical points of different nature. We show how various phases leave clear signatures 
to momentum distributions of the atoms which can be observed after time of flight expansion. 
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Recent advances in the experiments of ultracold Fermi 
gases have shown great potential for elucidating long- 
standing problems in many different fields of physics re- 
lated to strongly correlated Fermions. For instance, in 
recent experiments P, 0, H, 0, [B| spin-density imbalanced, 
or polarized, Fermi gases were considered. Such systems 
make it possible to study pairing with mismatched Fermi 
surfaces, potentially leading to non-standard phases such 
as that appearing in FFLO-states @, 0] or BP-states 
(Sarma-states). These possibilities have been consid- 
ered extensively in condensed-matter, nuclear, and high- 
energy physics The experiments in trapped gases 
have shown clear evidence of the separation of the gas 
into a BCS core region and a normal state shell around it, 
i.e. phase separation. Although an FFLO-type state has 
been predicted to appear in these systems as well OH EH, 
it is likely to be difficult to observe since it appears in the 
edges of the trap except for large polarizations. 

Recent experiments [13 . [l3| on Fermi gases confined in 
optical lattices have already demonstrated the potential 
of these systems for a multitude of studies of new phases, 
dimensionality effects, and dynamics. In this letter, we 
calculate the phase diagram for an attractively interact- 
ing Fermi gas in an optical lattice at zero and finite tem- 
peratures. In particular, we consider the possibility of 
the single mode FFLO-phase, where the order parame- 
ter is space-dependent, and of the BCS-BP phase, where 
compared to the standard BCS phase the excess polar- 
ization is carried by additional Bogoliubov excitations. 
We investigate their competition with the phase sepa- 
ration (PS) of the gas into normal and BCS superfluid 
regions. Our results reveal that a typical phase diagram 
as a function of polarization and temperature is as shown 
in Fig. [TJ phase separation is expected for small polar- 
izations and temperatures, whereas at zero temperature 
the FFLO state appears after some critical polarization. 
At a finite temperature, the FFLO-phase competes with 
the BCS-BP phase. As the temperature increases the 
BCS-BP phase becomes energetically favorable, and as 
the temperature is increased even further the BCS-BP 
phase gives way to a normal polarized Fermi gas. 

Furthermore, the phase diagrams reveal a Lifshitz 
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Figure 1: (Color online) The phase diagram of the imbalanced 
Fermi gas in a lattice. Colors: BP = blue (at P = this is 
actually BCS), FFLO = yellow, PS = red, normal = white. 
The average filling is 0.2 atoms/lattice site in each component, 
J = 0.07E R , and U = -0.26.Efl, where E R = h 2 k 2 /2m is the 
recoil energy, k = 2n/\ and A = f030nm. Here Ao means 
the gap at T = and P = 0. The T c at P = is 41 nK. 



point which is surrounded by the normal, FFLO, and 
BCS-BP phases. The transitions around this point are 
of second order, but the FFLO phase breaks the transla- 
tional symmetry which is present in BCS-BP and normal 
phases. Interestingly, there is also a point surrounded 
by the PS domain and the FFLO and BCS-BP phases. 
Comparing this point to the tricritical point discussed in 
Refs. 14|, [l5|, the difference is that the part of the phase 
diagram around the tricritical point which was occupied 
by the normal state in Refs. [3, EH], is now occupied 
by the symmetry breaking FFLO phase. If extending 
to stronger interactions, the order of the transitions in 



question should be considered more carefully. 

In atomic gases, the mismatch of the Fermi surfaces 
is realized by a fixed atom number difference between 
the two components, i.e. fixed polarization. This is in 
contrast to many other physical systems (such as super- 
conductors or quark matter) where it is caused by dif- 
ferent chemical potentials, due to e.g. an external mag- 
netic field. In such systems it is known that, at zero 
temperature, the BCS state is realized for small chem- 
ical potential differences, and after a critical difference, 
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some exotic pairing state such as FFLO or BP state may 
appear, althou gh t he stability of these states is under a 
debate 1(| 17, lq |. If the polarization, rather than the 
chemical potential difference, is fixed, the BCS state can- 
not be the solution at zero temperature. Consequently, 
it is important to consider the phase separated state as 
an alternative to the FFLO and BP states. In fact, it 
turns out that at zero temperature the phase separation 
into an unpolarizcd BCS region and a polarized normal 
region gives a lower energy than the BP phase [HI ]. 

The phase diagrams including the possibility of FFLO- 
phase and phase separation have been calculated for 
atomic Fermi gases in free space [3, Also there, 

like in our results, the phase separation appears for low 
polarizations and the FFLO phase for higher ones. How- 
ever, a striking difference is that our calculations in opti- 
cal lattices give a considerable parameter window for the 
existence of the FFLO phase, whereas in free space the 
parameter window is extremely narrow. Recently, den- 
sity polarized gases in optical lattices were considered 
in Ref. but that work concentrates on the insulat- 
ing phases and does not resolve the competition between 
FFLO, BCS-BP, and phase separation, moreover the dis- 
cussion there is limited to zero temperature. Finite tem- 
perature phase diagrams for trapped imbalanced gases 

. These works do 



the Fermi function and the coefficients u and v are defined 
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have been considered in 
not consider the FFLO state, except [llj. There a FFLO- 
type phase appears via the Bogoliubov-de Gennes ansatz, 
not via a plane wave ansatz customary in translationally 
invariant or periodic systems, therefore the spatial be- 
haviour of the order parameter is different: for small and 
intermediate polarizations the oscillations of the order 
parameter appear only in the edges of the trap. 

We consider two distinguishable species of fermionic 
atoms, e.g. atoms in different hyperffne states, labelled 
by up and down pseudospin index, interacting via an 
attractive contact interaction and confined in an opti- 
cal lattice. The system is described by the Hamilto- 
nian H = £k(E CT 0k - ii a )c\ 
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h.c.) — jj, where 



2 J2i=i M 1 ~ cos(M))- Here q 
is the FFLO momentum and it represents the momen- 
tum of the Cooper pairs. In the single mode ansatz 
the spatial dependence of the order parameter is given 
by A(r) = Aexp(2zq • r) . The parameters U and 
Ji are defined as in [231 ]. Our calculations are in the 
intermediate and weak coupling regime (U/J < 6), 
where the single band approximation is valid. With 
q = the Hamiltonian is the standard BCS-BP Hamil- 
tonian. The quasiparticle energies of the system are 
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A 2 , where 
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fi a . We consider fixed numbers of particles 
and determine the chemical potentials by the number 

and 



equations 7V T = J2w u l q /(^k,q, 

^ = E k < q /(-£ k , q ,' 
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H < q /(£k, q , 
■£'k,q,+), where / is 



in the usual way, see for example 2J] . The grand poten- 
tial of the system is f2 = ^ k 



£l-k+q + -Ek, q , 
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^Ek [m(l + e- /3i5k '0'+) +ln(l + e' 3 ' Bk ^-)]. For any 
given point in the phase diagrams, we look for the stable 
phase by minimizing the free energy F — 17+ ^n-\+ ^[n^ 
with respect to q and A, where n-\ and are the filling 
factors of the different atoms. We define the different 
phases as follows: the region of the phase diagram where 
the energy is minimized by A = is normal phase, the 
region where A > and q > is FFLO. The regions 
where A > and q — are BP, and within the BP re- 
gion, the line with no polarization, i.e. n-\ =714, is the 
standard BCS state. All of these phases are obtained as 
stable solutions from the free energy given above and the 
BCS state can be reached smoothly from the BP state 
by reducing the polarization to zero. 

We also compare all the above states to a phase sepa- 
ration between the standard BCS state and the polarized 
normal state. We do this as in Ref. [lj|: with numbers Nf 
and -ZVj in the different components, we let Nbcs atoms 
from both components occupy the BCS state and put the 
rest of the atoms to the normal state. We let the normal 
state occupy a fraction x of the lattice, leaving the frac- 
tion 1—x for the BCS state. We then add the free energies 
of a BCS state with a density UBCs/iX ~ x ) m both com- 
ponents and a normal gas with densities (nij — rises) /% 
and (ni — rises) /% m the different components. There- 
after we vary both Nbcs and x to find the energy of 
the optimal phase separated state and compare it to the 
optimal free energy given by a single state. 

Fig. [2] shows the zero temperature phase diagram as 
a function of polarization P and the average lattice fill- 
ing n, for two different interaction strengths. Starting 
from zero, the critical polarization for the FFLO-normal 
boundary increases as function of the average filling. The 
maximum critical polarization is achieved for a certain 
filling factor that is the smaller the stronger the interac- 
tion strength is, and is different from half filling. Note 
that although half filling is predicted to give the maxi- 
mum gap [25j . in our case the value 0.5 of the average 
filling n, for high polarizations, is far from the actual half 
filling. In fact, for large n and P we would enter a regime 
outside the scope of our Hamiltonian restricted to the 
lowest band.[28| The phase separation - FFLO bound- 
ary shows interesting behaviour: both critical polariza- 
tions (for the PS-FFLO and FFLO-normal) are close to 
each other for small fillings, then grow when the density 
increases, however, the PS-FFLO critical polarization 
starts to decrease earlier and stronger than the FFLO- 
normal critical polarization. This implies that, for high 
densities, the polarization window for the FFLO phase 
can become large compared to the polarization window 
for the PS domain. For small densities the critical polar- 
ization tends to zero and the window for FFLO phase be- 
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Figure 2: Zero temperature phase diagram as a function of 
average filling n and polarization P for two different interac- 
tion strengths. The transition lines are with U = ~0.36Er 
(solid) and U = — 0.26-Er (dashed). We indicate the actual 
data points and the lines are meant as guides to the eye. We 
choose J = 0.07-Eh, Eh is the same as before. 



comes narrow, as expected in the long wavelength limit 
where the lattice dispersion can be approximated with 
the free space dispersion. 

The stability of the FFLO state in a lattice is due to 
the flatness and nesting of the majority and minority 
component Fermi surfaces. In the continuum limit the 
Fermi surfaces are spheres, and the FFLO ansatz effec- 
tively shifts them relative to each other, allowing match- 
ing of parts of the Fermi surfaces. For small lattice fillings 
n the well known result q ~ fcpf — fcpj. gives a good esti- 
mate for the q we find numerically. For higher fillings the 
surfaces deform from spherical towards octahedra. By 
shifting two octahedra relative to each other, four whole 
faces of the minor component octahedron can be made to 
match the major component Fermi surface. Such nest- 
ing makes the stability of FFLO intuitively understand- 
able. For fillings close to 0.2 where the Fermi surfaces 
can be approximated by octahedra, we find that the nu- 
merically obtained q can be well estimated by the shift 
needed along the diagonal of the octahedra required to 
make the surfaces match. For intermediate cases between 
spherical and octahedron, we made numerical optimiza- 
tion of the match of the corresponding ideal gas Fermi 
surfaces and found again good agreement with the values 
of q given by the full FFLO calculation. Typical values 
given for qd by our FFLO calculations are between and 
0.4, larger fillings giving larger values of q. 

Fig. [T] shows the finite temperature phase diagram. 
The FFLO state dominates over the phase separation 
after some finite critical polarization is reached, even 
at zero temperature. At high enough temperatures the 
phase separation gives way to the FFLO, BCS-BP, or 
normal phases. Note that the FFLO regime, both at 
zero and finite temperature, is reasonably large, not van- 
ishingly small as in many other systems, which makes 
the observation of this phase in optical lattices more 



feasible. We have calculated finite temperature phase 
diagrams for several values of the interaction strength 
U and the average density n. In general, larger values 
of these parameters produce higher critical polarizations 
and critical temperatures for the FFLO and BP states. 
The PS-FFLO boundary has a non-trivial dependence 
on the density, as already indicated by Figure 2. Over- 
all, the qualitative behaviour presented in Figure 1 is 
not sensitive for the variation of the parameters. For 
the parameters of Figure 1, A(P = 0,T = 0) ~ 0.17E F 
and T C (P = 0) ~ O.lEp. Including pair fluctuations 
[2r3 | would shift the normal state boundary down in tem- 
perature at higher temperatures but would not change 
the qualitative features such as the large FFLO win- 
dow found in lower temperatures. We have checked that 
including Hartree corrections does not change this re- 
sult either. At T = 0, the polarizations P for which 
phase separation (PS) and FFLO appear are, for three 
different interaction strengths U / (6 J) = U /W , the fol- 
lowing: phase separation dominates for P ~ [0 : 0.23] 
(U/W = 0.61, gap (= A(P = 0,T = 0)) ~ 0.17 Bp, 
P ~ [0 : 0.08] (U/W = 0.42, gap - 0.06 E F ), P ~ 
[0 : 0.05] (U/W = 0.37, gap - 0.04 E F ), and FFLO for 
P ~ [0.23 : 0.34], [0.08 : 0.13], [0.05 : 0.08], respectively, 
from where the relative size of the FFLO domain versus 
the total PS+FFLO domain is obtained to be 0.32, 0.38, 
0.38, respectively, i.e. it is essentially the same also for 
smaller interactions (even slightly grows). 

All the phases that we discussed here are directly ob- 
servable in optical lattices from the momentum distribu- 
tions that are obtained by imaging after time of flight 
expansion. Also other methods, such as observing the 
noise correlations [27}can be considered. Fig. [3] shows 
how the FFLO (at 15 nK) and BP (at 20 nK) phases 
are reflected in the momentum distribution. A normal 
state would be completely symmetric in the fc-space, in 
contrast to the FFLO, and the phase separated state can 
be assumed to produce a symmetric result as well (with 
a background trap, the actual form of phase separation 
is likely to be a BCS state in the middle of the trap and 
the normal gas on the edges). Furthermore, the phase 
separated state can be distinguished, for example, by re- 
moving the paired atoms from the system using RF-ficlds 
(detuning matched to the pairing gap to select the paired 
atoms only) and then observing the remaining cloud of 
atoms which is, after the removal, expected to have a 
non-monotonous density distribution if the initial state 
was a phase separated state. 

The background trapping potential that is present in 
most optical lattice experiments would further decrease 
the energy of the phase separated solution because of the 
higher density in the middle of the trap. The background 
trap must therefore be sufficiently slowly varying for the 
FFLO state to remain energetically favourable. An or- 
der of magnitude estimate for the maximum trapping 
frequency can be made by considering the energy differ- 
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Figure 3: (Color online) Difference in the momentum distri- 
butions, n-f —ni, integrated over the z-direction, of the two 
components with an average filling of 0.2 atoms/site in the BP 
(left) and FFLO (right) states. Darker colour means larger 
difference. The q is marked on the k x axis in the FFLO case. 



ence between the FFLO and phase separation solutions, 
per particle, SE which for typical parameters used here 
is about 10~ 3 Er. We compare this to the potential en- 
ergy of a harmonic trap muj 2 x 2 /2 at the distance x given 
by the size of the system, i.e. x — Nd/2, where N is 
the number of sites in one dimension and d is the lattice 
spacing. This gives the maximum energy that is possible 
to save by moving a particle from the edge to the middle; 
we also multiply this by the polarization to correspond to 
the particles that would likely be affected by the phase 
separation. Comparing these two energies, for the pa- 
rameters used here ( 6 Li, d = 0.5fim, N — 64), we obtain 
^max of the order of 100Hz, which is not unreasonable. 

In summary, we presented for the first time a finite 
temperature phase diagram for optical lattices where var- 
ious forms of superfluidity and phase separation are com- 
pared. Apart from giving the scales for critical polar- 
izations and temperatures, relevant for the experiments, 
we find a clear qualitative difference to similar works 
in other systems: the parameter window for the exis- 
tence of FFLO can be rather large, the phase is certainly 
not negligible or an edge effect. It turns out that the 
FFLO phase often appears in the parameter region where 
the lattice problem cannot be approximated by the long 
wavelength limit, which explains the pronounced differ- 
ence with the free space results, where the FFLO phase 
is expected only in a very narrow parameter region. An- 
other remarkable difference between optical lattices and 
other systems realizing imbalanced Fermi systems is that 
the signatures of the different phases should be readily 
observable. We show that the FFLO and BCS-BP states 



are clearly reflected in the momentum distribution [24 1 
of the gas, which can be directly observed after a time 
of flight expansion. Our study has been restricted to a 
single mode FFLO ansatz and it is likely that multi-q 
FFLO superfluids have a lower energy. Therefore, re- 
gions of our phase diagrams where FFLO is expected to 
appear may be larger and also contain several different 
types of FFLO phases. 
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